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ABSTRACT 

We study the large-scale velocity fields traced by galaxy clusters in numerical 
simulations of a box of side 960 Mpc, and compare them with available 
data on real clusters. In order to test the reliability of the simulations, which 
are based on an optimized version of the Zel'dovich approximation, we com- 
pare their cluster velocities with those of "exact" N-body simulations, and find 
a remarkable agreement between the two according to a variety of statistical 
tests. We analyse Cold Dark Matter (CDM) models with density parameter 
in the range 0.2 < fio < 1, both with and without the cosmological constant 
term to provide a flat geometry. We also simulate a Cold+Hot Dark Matter 
(CHDM) model, with 30% provided by the hot component. Comparison with 
real data is performed by applying tests based on the cumulative velocity fre- 
quency distribution (CVFD) and bulk flow statistics. For the CVFD, we use 
observational velocity data from different authors. By merging all available 
data in a combined sample, the CVFD test is able to exclude only open mod- 
els with rig < 0.4 and flat models with Slo — 0.2. However, the analysis of 
individual observational samples gives contradictory results; in particular the 
recent Giovanelli (1995) data, which imply much lower velocities (with signifi- 
cantly smaller uncertainties) than the other data, are actually only consistent 
with the open fig = 0.4 and the flat fio — 0.2 models, i.e. the models which are 
excluded by the other data. Furthermore, CVFD and bulk flow analyses of the 
Branchini et al. (1995) reconstructed velocity data disfavour precisely those 
models accepted on the grounds of Giovanelli's sample. Finally, we confirm 
that the Lauer & Postman (1994) bulk flow determination is an extremely 
rare event in the cosmological models we have analysed. 

Key words: Cosmology: theory ~ dark matter - galaxies: clustering, forma- 
tion - large-scale structure of Universe 
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1 INTRODUCTION 

The quality and quantity of galaxy peculiar velocity 
data derived from the Infra-Red TuUy-Fisher (IRTF) 
and Dn-cr scaling relations is progressively increasing, 
providing am ever more precise description of large-scale 
cosmic flows. Moreover, an independent technique to es- 
timate the distance of high redshift galaxies with very 
small errors, based on the shapes of type la supernova 
light curves, recently developed by Riess, Press & Kir- 
shner (1995), could allow a description of motions on 
even larger scales. Consequently the study of large-scale 
flows is potentially one of the most powerful tools with 
which to understand the structure of the universe on 
large 50 Mpc) scales (see reviews by Dekel 1994; 
Strauss & Willick 1995 and references therein). Further- 
more, the hypothesis of gravitational instability allows 
one to relate the peculiar velocity field directly to the 
entire mass distribution, under the plausible assump- 
tion that at such large scales galaxies share the same 
dynamics as the dark matter distribution. 

Clusters of galaxies provide a particularly efiicient 
and precise way to estimate peculiar motions, since 
the determination of many redshift independent dis- 
tances of cluster members leads to a substantial reduc- 
tion in the statistical errors compared to studies of iso- 
lated galaxies. Moreover, new techniques based on the 
Sunyaev Zel'dovich (1980) effect (Haenhelt & Tegmark 
1995) and on the brightest cluster galaxy distance in- 
dicator method (Lauer & Postman 1994) are beginning 
to provide independent estimates of cluster peculiar ve- 
locities whose reliability is expected to increase consid- 
erably in the next few years. 

A further important motivation for considering 
clusters is that they can sample the cosmic velocity field 
up to distances larger than those presently accessible by 
galaxies. On these larger scales, fluctuations are largely 
in the linear regime and are, therefore, more closely re- 
lated to the "initial conditions" from which large-scale 
structure developed. In addition, since large-scale ve- 
locities are mostly sensitive to long-wavelength density 
fluctuations, one expects that, for a fixed choice of the 
fundamental cosmological parameters (Hubble constant 
Ho, density parameter Qo and cosmological constant 
term Oa), the amplitude of large-scale motions is al- 
most independent of the composition of the dark mat- 
ter (DM) component (i.e., cold vs. hot DM). For these 
reasons, several authors (e.g., Bahcall, Cen & Gramann 
1994; Croft & Efstathiou 1994; Bahcall, Gramann & 
Cen 1994; Cen, Bahcall & Gramann 1994; Strauss et 
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al. 1995; Gramann et al. 1995) have recently addressed 
the issue of whether cluster motions can be used to 
constrain cosmological models. Various statistical tests 
have been used for this purpose, among which the ve- 
locity frequency distribution, bulk flows, velocity dis- 
persions, pairwise velocities and the velocity correlation 
function feature prominently. 

The present work is devoted to a quantitative com- 
parison between presently available observational data 
on cluster peculiar velocities and an extended set of nu- 
merical simulations of a number of different cosmologi- 
cal models. The simulations are based on an optimized 
version of the Zel'dovich approximation (Zel'dovich 
1970; Shandarin & Zel'dovich 1989; Coles, Melott & 
Shandarin 1993), the ability of which to simulate the 
large-scale distribution of galaxy clusters has been as- 
sessed by Borgani et al. (1995; see also Plionis et al. 
1995). We will show in the following that the quasi- 
linear description of the large-scale velocity field pro- 
vided by our simulation method reproduces to a good 
accuracy that obtained from N-body simulations. At 
the same time, the method is so cheap, computation- 
ally speaking, that it enables us to consider easily a 
large number of realizations of a large number of differ- 
ent scenarios. 

In this paper we consider 57o < 1 CDM models both 
with and without a cosmological constant term A im- 
posing the condition of spatial flatness: Qa = A/{3Hq). 
We also consider, as reference models, the Qq = 1 CDM 
and CHDM scenarios. Our main aim here is to investi- 
gate whether the cluster velocity field can constrain the 
values of Qq and SIa- (In fact, for a fixed spectrum nor- 
malization, the resulting motions essentially depend on 
the value of SIq: the larger the density parameter, the 
larger are the peculiar velocities.) 

Observational motivations for cosmological scenar- 
ios with fio < 1 have been considered by different au- 
thors (cf. Coles & Ellis 1994; see Ratra & Peebles 1995 
for a summary; see, however, Primack 1995, for a dif- 
ferent view about a low-Slo Universe). Due to the fact 
that "ordinary" inflation models can accommodate open 
universes only by fine tuning either the duration of in- 
fiation or the pre-infiationary conditions (e.g. Lucchin 
& Matarrese 1985; Ellis, Lyth & Mijic 1991), interest, 
in the past, has been limited only to models in which 
f^A 7^ (Efstathiou, Sutherland & Maddox 1990; Kof- 
man, Gnedin & Bahcall 1993). Only very recently has it 
been realized that, in the framework of the bubble nucle- 
ation model, it is possible to obtain an open geometry. 



regardless of the initial conditions, but requiring only an 
appropriate choice for the scalar field potential (Sasaki 
et al. 1993; Bucher, Goldhaber & Turok 1995); see Linde 
& Mezhlumian (1995), and references therein, for a re- 
view of this class of infiationary models. An extensive 
comparison of the predictions of open CDM models with 
a variety of observational data is given by Liddle et al. 
(1995) and Yamamoto & Bunn (1995); see also Ratra & 
Peebles (1994). 

The layout of this paper is as follows. In Section 2 
we present our approach for simulating the cluster ve- 
locity field and we discuss its reliability by comparing 
relevant results with those obtained using an N-body 
code. We also describe the statistical tests applied in 
the analysis and briefiy introduce the suite of cosmolog- 
ical models we consider. In Section 3 we compare the 
output for different initial spectra in terms of both the 
velocity frequency distribution and the velocity correla- 
tion function. Section 4 is devoted to a comparison with 
available observational data; in particular we study the 
bulk fiow and the cumulative velocity frequency distri- 
bution. In Section 5 we discuss our results and, finally, 
we draw our main conclusions in Section 6. 

2 THE SIMULATIONS 

2.1 The Zel'dovich approach 

Our simulations mimic gravitational dynamics through 
the Zel'dovich approximation (ZA; Zel'dovich 1970; 
Shandarin & Zel'dovich 1989). This method is based 
on the Eulerian-to-Lagrangian space mapping 

x(q,t) = q-b(t)V^(q), (1) 

where q and x are initial and final particle positions, 
respectively, b{t) is the fluctuation linear growth factor 
and the velocity potential ■(/'(q) is related to the ini- 
tial fluctuation fleld by the Poisson equation. The sim- 
ulation procedure has been described in detail in Bor- 
gani et al. (1995; hereafter Paper I), to which we refer 
the reader for more details. Here we simply recall that 
the initial power spectrum is convolved with a Gaus- 
sian window, P{k) P(k)e~'' , in order to suppress 
shell-crossing at small scales (Coles, Melott & Shan- 
darin 1993; Melott, Pellman & Shandarin 1994). The 
filtering radius Rf is chosen in such a way that the av- 
erage number of streams at each Eulerian point (Kof- 
man et al. 1994) is Ns — 1.1. After assigning the lin- 
ear potential on a grid, particles are moved from their 
initial grid positions according to eq.(^. The density 
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and potential fields are then reassigned on the grid and 
clusters are selected as the highest local density max- 
ima, so as to reproduce the observed AbcU/ACO cluster 
number density, rid — 1-8 x 10^'''( ft^^Mpc)"^, which 
corresponds to dd ^ 38/i~^Mpc for the mean cluster 
separation. The parameter h is the Hubble constant in 
units of 100 km s"^ Mpc"^ 

Velocities at the grid positions are evaluated by fol- 
lowing two different prescriptions. 

(a) Defining the particle velocities as v oc — V'(/'(ci) 
and interpolating their mass and momentum on the grid 
with the TSC scheme (Hockney & Eastwood 1981). The 
velocity at the grid point is then defined as the ratio 
between the local values of momentum and mass (e.g., 
Kofman et al. 1994). This linear prescription, that we 
will name LIN in the following, amounts to the assump- 
tion that particles at the final positions do not still feel 
any tidal force and move according to the initial {linear) 
gravitational potential. 

(b) Computing the gradient of the ZA-evolved poten- 
tial, > which is connected by the Poisson equation 
to the density fluctuations, , traced by the final par- 
ticle distribution, V'^ip^A ~ -iT^Gpa^S^j^. Accordingly, 
V oc — Vi/'zyi. In this quasi-linear prescription, named 
ZEL in the following, one assumes that, although the 
density field has undergone a substantial non-linear evo- 
lution, it is still connected to the velocity by a linear 
relationship. 

Non-linear effects in the velocity field, such as in- 
fall and merging, are not accounted for in this simula- 
tion method. These effects, however, are only expected 
to dominate on scales ^ 10-20 ft" ^Mpc, where the ob- 
servational data on cluster peculiar velocities are in any 
case rather unreliable. 



2.2 Comparison with N— body velocities 

In Paper I we noted that, while the ZA is very accurate 
at locating clusters in the correct positions for models 
with (Tg^ 1 [ug is the rms fluctuation amplitude within a 
top-hat sphere of radius 8/i^^Mpc radius], non-linear 
gravitational effects on the cluster mass scale degrade 
this precision when crg^ 1. However, this is not expected 
to represent a limitation for the present analysis. In fact, 
in order to study the cluster velocity field it is not neces- 
sary to know cluster positions with any great precision. 
What we need from simulations is only a population of 
objects tracing the large-scale flows having the same 



sampling density and same selection biases as real clus- 
ters do. 

Before presenting the results of our analyses, it is 
important to assess the reliability of our approach to the 
simulation of the cluster velocity field. We therefore be- 
gin by comparing results of our approach with those ob- 
tained from the same N-body simulation we previously 
used in Paper I for the c:omparison of cluster positions. 
The initial spectrum corresponds to the Cold-I-Hot DM 
(CHDM) model with 30% of hot component (see next 
section). The box size is 320/i~^Mpc (with h = 0.5) 
and the number of grid points and particles is 128^. We 
do not distinguish between hot and cold particles, since 
any effect of residual free-streaming should be negligi- 
ble at the smallest scales 2.5 /i~^Mpc allowed by our 
resolution. 

In the N-body simulation clusters, are identified as 
local maxima on the grid, following the same method 
outlined in the previous subsection. The velocity at the 
grid point positions is defined by following prescription 
(a), previously described, with the difference that non- 
linear (final) particle velocities are now used and that 
density and momentum are smoothed with a top-hat 
filter of width 10 /i~^Mpc to create continuous fields and 
to ensure that the density is non-zero at each grid point. 

For comparison we run a ZA simulation having the 
same initial spectrum and random phase assignment 
and assuming the same box size and resolution. We show 
results at two different evolutionary stages, correspond- 
ing to ag = 0.67 and erg = 1. Even though the two-year 
COBE data (Bennett et al. 1994) are roughly consis- 
tent only with the first epoch, we consider the second, 
higher, normalization as well, in order to assess the reli- 
ability of the ZA approach in a more evolved situation, 
when the cluster mass scale starts to become non-linear. 
The filtering radii, Rf, applied in the ZA simulation are 
Rf = 2.3 and Rf = 4.5/i"^Mpc for ag = 0.67 and 
erg = 1, respectively. In the following we will compare 
the velocity fields traced by the clusters identified in 
the N-body simulation and in the ZA simulation, by 
adopting the two definitions (LIN and ZEL) of velocity 
described above. 

2. 2. 1 General features of the velocity fields 

In Figure 1 we compare the velocity fields traced by 
clusters within a slice 60/i~^Mpc thick for the PM and 
ZA simulations. The fields are recovered by smoothing 
the velocities associated to the PM and ZA cluster po- 
sitions with a Gaussian window of radius 20 /i~^Mpc. 
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The resulting velocity field is then reassigned on 16'^ 
grid points. All clusters are equally weighted, indepen- 
dently of their masses. At the first stage {as = 0.67), 
there is a strong correspondence between the cluster ve- 
locity fields obtained by different methods: all the main 
features are well reproduced at the correct locations. At 
as ~ 1, although the principal stream lines are still well 
delineated, there is a slight tendency for both ZEL and 
LIN to overestimate the velocity amplitude. 

It is quite interesting to follow the evolution of the 
non-linear feature at the centre of the panels displayed, 
where some degree of infall (convergence of stream lines) 
takes place. At erg = 0.67 this feature is quite well re- 
produced by the LIN and ZEL velocities. However, at 
as ~ 1 the higher degree of infall displayed by the N- 
body velocities is not reproduced by the LIN and ZEL 
fields. As the degree of non-linearity increases, shell- 
crossing becomes more and more important: stream 
lines consequently cross each other, instead of simply 
converging towards potential minima. 

2.2.2 Statistical tests 

In order to obtain a more quantitative assessment of 
the reliability of our approach, we analyze PM, LIN and 
ZEL velocity fields with the same statistical tests to be 
applied in the analysis of the data we perform later. 
They are the 3D velocity frequency distribution (VFD), 
the velocity correlation function and the bulk fiow (as 
a function of scale). Figure 2 summarizes the results at 
the two evolutionary stages. 

The VFD P{v) is defined as the fraction of cluster 
velocities in a given range. The results of P{v) for the 
3D cluster velocity field are shown in the upper panels. 
We used a Kolmogorov-Smirnov (KS) test to compute 
the significance level for the null hypothesis that the 
different VFDs are drawn from the same distribution. 
Comparing LIN to ZEL and PM, we find significant de- 
partures at both output times: the distribution peak for 
LIN is always shifted to lower velocities (~ 300 ver- 
sus ~ 450 — 500kms~^ at as = 0.67 and ~ 500 versus 
~ 600 — 650kms~^ at as — 1). On the other hand, the 
results for PM and ZEL are in good agreement: the KS 
test indicates high probabilities both for both as = 0.67 
and as = 1. 

Note that the VFD statistics do not take into ac- 
count the spatial distribution of clusters; at the more 
evolved epoch the agreement of cluster positions be- 
tween PM and ZEL is not as good (see the discussion 
in Paper I). In order to take into account the cluster 



positions in the statistics of the velocity field, we will 
use the velocity correlation tensor, defined as 

■9.,{r) = {v,{r,)v,{r2)). (2) 

Here, ri and r2 = ri -I- r are the position vectors of 
two objects having peculiar velocities v(ri) and v(r2), 
respectively (Gorski 1988; Gorski et al. 1989). In gen- 
eral what one estimates is the scalar velocity correlation 
function £,v{r) which is defined as the trace of the ve- 
locity correlation tensor, ^„ = I *I'ii|^^'^- According 
to its definition, ^(u) is a measure of the coherence of 
the velocity field: if at some scale a coherent bulk fiow 
dominates the velocity field, will be positive; infall 
or outflow velocities result in a negative contribution to 

The velocity correlation function (r) is plotted in 
the central panels of Figure 2. Error bars, shown only 
for the PM results, are la uncertainties computed using 
the bootstrap technique. The agreement between differ- 
ent determinations is good at both epochs: only at very 
small distances (r^ 20/i~^Mpc), some differences (al- 
ways smaller than 2a) are present. 

The last test we consider is the bulk velocity Vbuik, 
defined as the centre-of-mass velocity of a specified re- 
gion and given by the integral of the peculiar velocities 
v(x) over a selected volume specified by a selection func- 
tion (;/>(x): 

Vbuik(r)= / 0(x)v(x)dx. (3) 

"'0 

In our case the velocity field is traced by galaxy clusters 
and eq.(^ becomes 

Vbuik = ^ WiV,/ ^ Wi, (4) 

i i 

where i refers to different clusters inside the considered 
region. The weights Wi account for cluster masses, in- 
homogeneous sampling and uncertainties in the cluster 
positions. The choice of the weighting scheme is par- 
ticularly important when one deals with observational 
biases, and several different schemes have been pro- 
posed in literature. For the purpose of this analysis, we 
are studying a perfectly sampled cluster distribution, 
for which the cluster masses and peculiar velocities are 
known, so the natural weighting scheme to implement 
is simple number weighting (i.e. Wi = 1). 

Results for the bulk velocity statistics are shown in 
the lower panels of Figure 2. Also in this case we plot 
the error bars only for the N-body results: they are ob- 
tained as standard deviation of bulk ffows computed for 
10,000 different observers. Note that the rather small 
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box size is expected to result in a significant underes- 
timate of Vbuik at scales ^ 70-100 /i~^Mpc. However, 
here we are interested only in a relative comparison of 
the different simulation methods rather than estimating 
the absolute value of the different model bulk flows. At 
(Tg = 0.67, ZEL simulations provide a reliable estimate 
of Vbuik on scales larger than 20 /i~^Mpc. The situation 
is even better at erg = 1; here the agreement between 
PM and ZEL is well inside the la range at all scales. 
On the other hand, the LIN simulation tends to under- 
estimate the bulk flow at small scales when ag, = 1- 

We conclude from this analysis that the quasi- 
linear definition ZEL gives a good representation of the 
velocity field traced by clusters over the whole scale 
range, even at the stage erg = 1, when substantial non- 
linearity appears on the cluster mass scale. This is par- 
ticularly relevant for those statistics incorporating both 
position and velocity information, such as the velocity 
correlation function and bulk flow. We will therefore use 
the ZEL velocity definition in our cluster simulations 
from now on. 

2.3 Models for the power spectrum 

As we have explained, our aim is to use the cluster ve- 
locity field to constrain the value of the density param- 
eter Q,(i, both with and without a cosmological constant 
term, JIa, assuming the context of a CDM scenario. We 
considered values of the density parameters correspond- 
ing to Q.0 — 0.2, 0.3, 0.4, 0.6, 0.8. In addition, we also 
simulated two f^o = 1 models, namely the standard 
CDM model (SCDM) and a Cold+Hot DM (CHDM) 
model with 30% of hot component. We therefore end 
up with 12 models for the forthcoming analysis. They 
are listed in Table |l|, in which OCDM stands for the 
open models and ACDM for the models with a non- 
zero cosmological constant term. 

For the CDM power spectra we assumed the expres- 
sion provided by Efstathiou, Bond & White (1992), with 
the shape parameter F — QQhexp{—QB — /f^o), cor- 
rected according to the prescription of Peacock & Dodds 
(1994) to account for the presence of baryons (Sugiyama 
1995). We take the density parameter in baryons to be 
Q,B ~ 0.013/i^, as suggested by standard nucleosynthe- 
sis (e.g. Reeves 1994). The transfer function provided 
by Holtzman (1989) is used for the CHDM spectrum. 
For each pair of (f^o, ^a) values, we choose the Hubble 
parameter h in such a way that the age of the universe 



1 X y^Qox-^ + (1 — fio — flA)x'^ + SIa 



(see, e.g., Peebles 1993) is to — 12 Gyr for all the models, 
except for SCDM and CHDM, for which we assume h = 
0.5, giving to — ISGyrs. Although to — 12 Gyrs seems 
to be slightly too small in comparison to most current 
estimates of the age of the universe (e.g. Chaboyer et al. 
1995; arguments in favour of io — 11 Gyrs from globular 
cluster ages have also been advanced, e.g., Shi 1995), we 
adopt this value of to for the low-fio models so as to 
obtain almost the same age as with the usual choice of 
h = 0.5 for the Slo = 1 models. 

To normalize the two f2o = 1 models to COBE, we 
use Qrma-ps = 20fiK (Gorski et al. 1994). In the case 
of the OCDM models, we resort to the normalization 
procedure outlined by Gorski et al. (1995). The pro- 
cedure is straightforward to implement, since we select 
models with the same criteria as Gorski et al. (1995) did 
in their Table 1. We take for ag the central value in their 
reported la range, which is associated to uncertainties 
in the Qrms-ps estimates. 

For the ACDM models, Efstathiou et al. (1992) pro- 
posed the relation 

Q™.-ps= (g^)"(f)Vonr^iJ°% (6) 

(To — 2.735 K) which connects the purely Sachs-Wolfe 
(SW) quadrupole to the amplitude B of the primordial 
power spectrum, P{k) = Bk. However, as, for exam- 
ple Stompor, Gorski & Banday (1995; hereafter SGB) 
have pointed out, signiflcant additional contributions to 
the radiation power spectrum come from the integrated 
SW effect induced by the cosmological constant. Unfor- 
tunately, these authors reported, in their Table 2, the 
as normalizations only for models with two fixed val- 
ues of h, namely h — 0.5 and h = 0.8. In order to fix 
the normalization for the models corresponding to dif- 
ferent choices of h, we proceed as follows. SGB noted 
that the integrated SW effect is negligible for f2o?^ 0.5. 
Therefore, we use for Qq > 0.5 the pure SW normaliza- 
tion provided by eq.(^. On the other hand, from Figure 
3 of SGB, it turns out that for our Slo < 0.5 models 
the as normalization essentially depends on the value 
of ^loh^. Therefore, we compute the floh^ value for each 
model and use this plot to normalize our spectra when 
Qo < 0.5. We expect the accuracy of this procedure to 
be within 10%-15%, of the same order of the statistical 
and systematic deviations estimated by SGB in their 
analysis. 

The values of erg for all the considered models are 
shown in Table |^, where we also report the filtering radii 
Rf used to suppress shell-crossing in the ZA (see Sec- 
tion 2.1). 
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Table 1. The model parameters. Column 2: the density pa- 
rameter Qq; Column 3: the cosmological constant term Q/^; 
Column 4: the Hubble parameter h; Column 5: the linear 

rms fluctuation amplitude at 8/i~^Mpc as; Column 6: the 
filtering radius Rf (in h~^Mpc). 



Model 


no 


Qa 


h 


0-8 




SCDM 


1.0 


0.0 


0.50 


1.36 


6.0 


CHDM 


1.0 


0.0 


0.50 


0.78 


3.3 


OCDM02 


0.2 


0.0 


0.70 


0.31 


0.5 


OCDM03 


0.3 


0.0 


0.65 


0.50 


1.6 


OCDM04 


0.4 


0.0 


0.65 


0.75 


3.1 


OCDM06 


0.6 


0.0 


0.60 


1.10 


4.9 


OCDM08 


0.8 


0.0 


0.55 


1.30 


5.8 


ACDM02 


0.2 


0.8 


0.87 


1.20 


5.8 


ACDM03 


0.3 


0.7 


0.78 


1.44 


6.9 


ACDM04 


0.4 


0.6 


0.72 


1.-54 


7.2 


ACDM06 


0.6 


0.4 


0.64 


1.17 


6.6 


ACDAIO.S 


o.s 


0.2 




l.i:-! 


6.3 



For each model we ran a single realization within 
a box of 960/i~^Mpc on a side, using 256^ grid points 
and particles. With the adopted cluster mean separa- 
tion, we end up with about 16,000 objects in each sim- 
ulation box. The very large size of the box ensures that: 
(a) no appreciable box-to-box variance is present in 
the final results (we verified this by running two real- 
izations for ACDM02 and CHDM models); (b) fluctu- 
ation modes with wavelengths larger than the box size 
make a negligible contribution to bulk flows and veloc- 
ity correlations. Indeed, our analysis is confined within 
180/i~^Mpc, which is roughly the scale where the re- 
constructed bulk flow traced by Abell/ACO clusters is 
reliable (Branchini, Plionis & Sciama 1995; BPS here- 
after). 

3 COMPARISON BETWEEN MODELS 

In this section we consider the 3D velocity frequency dis- 
tribution and the velocity correlation function, with the 
aim of quantifying intrinsic differences between models. 
Any comparison with observational data and how ob- 
servational uncertainties are taken into account in sim- 
ulation analysis will be discussed in the next section. 

3.1 The velocity frequency distribution 

As a flrst test, we consider the VFD P{v) for the 3D 
cluster velocity field. We compute the VFD for the 
whole box for each of the simulations. In Figure 3 we 



plot the results for SCDM and CHDM models (top left 
panel), the OCDM models (top right panel) and ACDM 
models (bottom left panel). As expected, models with 
low Qo generate smaller velocities, independently of the 
presence of the cosmological constant. However, the f^A 
term tends to increase the average peculiar velocity and 
is also effective in concealing differences in the VFDs 
among various fio choices. Consequently, ACDM models 
have VFDs that are similar to each other for Qo > 0.4, 
while the analogous OCDM models show much better 
separated VFD curves. The model with the largest ve- 
locities is SCDM while the probability of having very 
large velocity in CHDM is smaller than in several of the 
low density models. 

Due to their small size, we prefer not to plot the 
error bars, estimated as scatter between the two realiza- 
tions of CHDM and ACDM02 models. Their smalhicss 
is a consequence of the large dimension of the simulation 
box and confirms that the box size we adopted is large 
enough to avoid any significant effect of cosmic variance. 
This is especially true for CHDM, which has a smaller 
relative amount of large-scale power and, consistently, 
displays a smaller dispersion than ACDM02. Therefore, 
we expect that uncertainties in the analysis of realistic 
simulated samples (see next section) will be associated 
with the limited number of included clusters and to the 
observer-to-observer scatter, rather than to the cosmic 
variance arising from the finite box size. 

Considering N body simulations of different cos- 
mological scenarios (SCDM, ACDM03 and Oo = 0.3 
primeval baryonic isocurvature models), Bahcall, Cen 
& Gramann (1994) showed that the shape of the cluster 
velocity distribution is well matched by a Maxwellian 
distribution P{v) oc cyip{—v'^ /2a^), where ct„ is the 
velocity dispersion. This result is expected when a Gaus- 
sian distribution for the initial density fluctuations is as- 
sumed and effects of non-linear gravitational clustering 
are negligible at determining cluster velocities. Using 
the value of cr„ obtained from the VFDs of our simula- 
tions, we found that a Maxwellian distribution provides 
a good representation for all considered models. In Fig- 
ure 3 the fitting curves are plotted (dashed lines) for 
clarity only for the two Oo = 1 models: they are almost 
indistinguishable from the simulated data and closely 
resemble the result of Bahcall et al. (1994a), which is 
based on N-body experiments. 
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3.2 The velocity correlation function 

Cen, Bahcall & Gramann (1994) showed that the veloc- 
ity correlation function £,v{r), especially at small separar 
tion, is strongly dependent on cosmological models, and 
therefore represents a potentially useful tool for con- 
straining them. However, with the observational data 
currently available it is very difficult to obtain a de- 
termination of the cluster velocity correlation function, 
which is sufficiently reliable to be effectively compared 
it with simulated data, although Croft & Efstathiou 
(1995) have recently succeeded in applying a maximum 
likelihood analysis to the velocity correlations estimated 
from two sets of observational data, finding marginal ev- 
idence for more power than the standard CDM model. 

In Figure 4 we plot ^u(r) for the considered mod- 
els at scales r > 20/i~^Mpc, where the ZA description 
of the velocity correlations has been shown to be reli- 
able. Consistent with the results from the VFD analy- 
sis, OCDM models are characterized by a progressively 
larger ^«(r) as Jlo increases. A clear discrimination be- 
tween models is possible up to scales r ~ 60/i~^Mpc, 
while above 100/i~^Mpc all the models seem to con- 
verge, with the exception of OCDM02. Among the 
ACDM models, the only one with a significantly differ- 
ent behaviour is ACDM02: all the other spectra provide 
almost indistinguishable velocity correlations. 

We note also that the velocity correlation functions 
of different reaUzations for CHDM and ACDM02 mod- 
els are extremely similar, thus indicating that cosmic 
variance is negligible in such big simulation boxes. 

4 COMPARING MODELS AND 
OBSERVATIONS 

In this section we attempt to constrain the different cos- 
mological models by comparing predictions of properties 
of the cluster peculiar velocity field with relevant avail- 
able data. To this end, we compute both the bulk flow 
and the cumulative velocity frequency distribution. We 
prefer not to consider velocity correlations, since esti- 
mating in a meaningful way would require a large, 
statistically complete sample of accurate cluster pecu- 
liar velocities, which is not available at present. 

4.1 Cluster velocity data 

Unlike the simulation case, estimating the bulk flow ex- 
ecuted by real clusters is a non-trivial task, since the 
available cluster IRTF and Dn-a peculiar velocity data 



are sparse and inhomogeneous. Furthermore, only the 
line-of-sight component of the peculiar velocity is ob- 
served. These limitations can signiflcantly affect the de- 
termination of the observational bulk velocity (Kaiser 
1988; Regos & Szalay 1989) and they have to be taken 
properly into account if one wants to compare observa- 
tional and model bulk flows. Therefore, in the following, 
we will use the IRTF and D„ a cluster peculiar veloci- 
ties to compute only the cunmlative velocity frequency 
distribution (CVFD), P{> v), i.e. the fraction of clus- 
ters with a velocity greater than v (the VFD is its differ- 
ential version) and compare it with model predictions. 

A possible way to avoid the above-mentioned prob- 
lems in the context of the bulk flow test is to reconstruct 
the full 3D cluster peculiar velocity fleld from their ob- 
served redshift space distribution [see Dekel (1994) and 
Strauss & Willick (1995) for reviews of the reconstruc- 
tion methods available]. In the following analysis we 
will consider the bulk velocity derived from the recon- 
structed Abell/ACO cluster density field of Branchini 
& Plionis (1995) and Branchini et al. (1995), hereafter 
BPS. These authors found a good consistency between 
their bulk flow and their derived 3D smoothed density 
fleld, and independent measurements contained within 
a depth of 8000 km s~^ in the Mark III galaxy sam- 
ple (cf Dekel 1994; Hudson et al. 1995). We will also 
use the BPS velocities to constrain model predictions of 
P{>v). 

4.1.1 Observational velocities 

We construct our composite observational catalogue 
from cluster peculiar velocity determinations by Aaron- 
son et al. (1986), Mould et al. (1991,1993), Han & Mould 
(1992) and Mathewson, Ford & Buchhorn (1992), which 
are all based on IRTF measurements. We add the fur- 
ther data reported in Table 7 by Mould et al. (1991), 
which contains peculiar velocities obtained both with 
the IRTF method (Aaronson et al. 1989) and with the 
D„-a technique (Lucey & Carter 1988; Faber et al. 
1989). Finally, we include cluster velocities obtained by 
Hudson (1994), both with D„-a and IRTF techniques, 
and the recent IRTF data from Giovanelli (1995; see 
also Giovanelli et al. 1995). Of course, it is possible that 
more than one velocity determination is available for 
any given cluster: in such cases we have chosen to in- 
clude in the catalogue only the determination with the 
smallest error. The result of this selection is an interim 
sample of 65 objects, with a very inhomogeneous spa- 
tial distribution. Some of these clusters show very large 
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Table 2. The composite catalogue of observed cluster pe- 
culiar velocities. Column 1: the source of the data; Column 
2: the method used; Column 3: the number of cluster pecu- 
liar velocities available in the original paper; Column 4: the 
number of objects actually used in the composite catalogue. 



Source 


Method 


Available 


Used 


Aaronson et al. (1986) 


IRTF 


10 





Mould et al. (1991) " 


IRTF 


9 


6 


Mould et al. (1991) 


IRTF 


16 


1 


Mould et al. (1991) " 


D„-a 


16 


1 


Han & Mould (1992) 


IRTF 


21 


6 


Mathewson et al. (1992) 


IRTF 


24 


10 


Mould et al. (1993) 


IRTF 


8 


8 


Hudson (1994) 


IRTF 


21 


1 


Hudson (1994) 




17 


6 


Giovanelli (1995) 


IRTF 


23 


23 


TOTAL 






62 



" from Table 6. 

from Table 7 (Aaronson et al. 1989). 
= from Table 7 (Lucey & Carter 1988; Faber et al. 1989). 



errors and we prefer to exclude them from our analy- 
sis to obtain a more robust result. We therefore decided 
to consider only objects with peculiar velocity uncer- 
tainties smaller than 850 km s~^, excluding 2 objects, 
namely 2159-32 from Mathewson et al. (1992) and Her- 
cules from Aaronson et al. (1986). We have verified that 
our conclusions do not change if we exclude the clusters 
with velocity errors larger than 450 km s~^, so as to 
end up with a smaller catalogue of 43 objects. Finally, 
we also reject the data of Cen45, whose determination 
is uncertain and controversial, with different authors re- 
porting greatly discrepant results [i.e. 203 ±383 km s~^ 
from Han & Mould (1983), while 1663±336 km s"^ from 
Mould et al. (1991); note that both results are obtained 
with the IRTF method]. At the end, the composite sam- 
ple contains 62 objects, whose sources are detailed in 
Table |. 

Due to the two different methods used (IRTF and 
Dn-a relations) and to the absence of an inner cali- 
bration of the data coming from different sources, our 
composite sample may be affected by some basically un- 
controllable biases. For these reasons, after discussing 
results coming from the composite sample, we will also 
present results obtained using data coming from unique 
sources. 



4.1.2 BPS reconstructed velocities 

Branchini et al. (1995) adopt a two step procedure to 
recover the real-space positions and peculiar velocities 
of Abell/ACO clusters, within cz < 25,000 km 
starting from their redshift space positions. 

• As a first step, BPS generate, via Monte Carlo tech- 
niques, several different realizations of a synthetic clus- 
ter population, spatially correlated with the observed 
clusters, so as to correct for selection effects and prop- 
erly homogenize the Abell and AGO catalogues. The 
outcome is an all-sky volume limited distribution of syn- 
thetic + real clusters in redshift space. 

• The second step is the application of an iterative 
technique, similar to that proposed by Strauss & Davis 
(1988) and Yahil et al. (1991), which assumes linear 
theory and linear biasing, to reconstruct the true cluster 
positions from their observed redshifts. 

The velocity of each cluster is, thus, essentially ob- 
tained from its measured dipole evaluated at its recon- 
structed real space position. Taking into account the un- 
certainties due to the intrinsic error of their method as 
well as the uncertainties related to modelling the obser- 
vational parameters (galactic absorption, Abell/ACO 
homogenization procedure, etc.), BPS estimated their 
mean cluster velocity error to be ~ 170 km s^^ with a 
dispersion of about 90 km s~^ around this value. 

We note that the BPS reconstruction method is 
based on two assumptions, namely that: (a) linear grav- 
itational instability holds and (b) the simple linear bi- 
asing relation, Sn/n — b(5p/p), connects fluctuations in 
cluster number counts and DM density, with a position- 
independent biasing parameter b. While the former as- 
sumption is expected to be satisfied at the large scales 
relevant in our analysis, the latter is probably an over- 
simplification. Although its effect, in the fight of the 
success of the reconstruction method (cf. BPS), is likely 
not to be dramatic it could nevertheless introduce un- 
certainties into the individual cluster velocities which 
are not included in the previous quoted errors. 

We wiU use the 280 clusters within cz < 20, 000 km 
s~^, since BPS consider their cluster velocity field to be 
reliable out to to this depth. Beyond cz = 20, 000 km 
s~^ the redshift selection function drops exponentially 
while for cz > 25, 000 km s~^ they have used the rather 
crude assumption of a homogeneous universe in their 
reconstruction algorithm. 
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4.2 The Bulk Flow Analysis 

In order to perform a consistent comparison between 
simulations and real data, one would desire in princi- 
ple to reproduce the observational set-up as closely as 
possible. We place observers at grid points and define 
their local density contrast Sohs and bulk velocity aver- 
aged over a top-hat sphere of radius R — 7.5/i~^Mpc 
centred on them. The local value of S^hs is estimated 
by convolving the Fourier transform of the fluctuation 
density field with the Fourier transform of the top-hat 
window, W{kR) = 3(sinfcii- kRcoskR) / {kRf . The 
observer's velocity is defined using the same procedure 
adopted for the cluster velocity. We estimate Jobs and 
the peculiar velocity for 80,000 observers chosen at ran- 
dom grid positions in each simulation. We then select 
those 'observers' that have the same characteristics as 
the Local Group, the number of which is also considered 
as a test for DM models (Gorski et al. 1989; Tormen et 
al. 1993; Strauss et al. 1995; Moscardini et al. 1995; Tini 
Brunozzi et al. 1995). These characteristics are: 

(1) peculiar velocity Vlg = 627±44kms~^ (error cor- 
responding to 2(T uncertainties; Kogut et al. 1993) for 
a top-hat sphere of radius R — 7.5/i~^Mpc centred on 
the observer; 

(2) density contrast within the same sphere in the 
range —0.2 < 5lg < 1- 

We decided to neglect in this analysis a further require- 
ment concerning the quietness of the local flow, which 
implies a small value for the local shear. In fact, we ex- 
pect that any shear is poorly represented by the velocity 
description provided by our simulations. In any case, as 
shown by Moscardini et al. (1995), the local shear con- 
straint is not a particularly restrictive one. 

In Column 2 of Table ^ we report the fractions Tlg 
of observers satisfying the LG requirements. The results 
for different models are quite similar, with the highest 
probability occurring for OCDM06 and ACDM03 mod- 
els. The only exceptions are represented by the open 
models with f2o < 0.3, where none of the chosen grid 
points has the desired characteristics. This problem for 
the open models with very low fio and small ag was al- 
ready noted in Tormen et al. (1993). Due to the impos- 
sibility of deflning LG observers in these cases, we prefer 
to consider in the following analysis 2000 observers ran- 
domly chosen in the list of 80,000, dropping completely 
the LG constraints. We checked that, at least for those 
models having a large enough number of LG-like ob- 



servers, this choice does not affect the results signifi- 
cantly. 

A self-consistent comparison with the BPS bulk ve- 
locity would, strictly speaking, require us to apply their 
method also on the simulated cluster distribution, start- 
ing from their redshift -space positions. We have per- 
formed such a test on a subset of CHDM and ACDM 
observers and we found a very good agreement between 
the profiles of the "true" and reconstructed bulk flows. 
On the other hand, their amplitude can be matched a 
posteriori, by an appropriate choice of the (observer- 
averaged) cluster biasing parameter h (we note that 
such definition of b does not necessarily coincide with 
that based on the ratio of rms fluctuations in cluster 
counts and DM density; but see Kolokotronis et al. 
1995). Having assessed the agreement between "true" 
and reconstructued bulk flows, we decided to compare 
the BPS data to the original ZEL cluster velocities, since 
applying the reconstruction procedure for all the 2000 
observers selected in each model would require a very 
time-consuming computation. 

The cluster bulk velocity for the simulations was 
computed as follows: for each selected observer we mea- 
sure the bulk velocity in spheres of increasing radius 
using eq.(^) with Wi — 1 (we verified that very similar 
results are obtained using a mass weighing scheme, i.e. 
with Wi oc nii, where rrii is the cluster mass). Due to 
the large intercluster separation, we started measuring 
the bulk velocity within spheres of radius 50 h~^Mpc to 
have reliable estimates. A limiting depth of 180 /i~^Mpc 
has been chosen to allow a meaningful comparison with 
the BPS data. 

The scale-dependence of the bulk fiow averaged 
over the 2000 observers is plotted in Figure 5 for all the 
cosmological models we have considered. The shaded re- 
gion represents the la strip allowed by the BPS data. A 
simple visual analysis reveals that the only models we 
can rule out at a level larger than la are the OCDM 
ones with Qo < 0.4. Flat models with non-vanishing 
cosmological constant, as well as CHDM and SCDM 
are in good agreement with the data on scales larger 
than ~ 100/!.~^Mpc. On smaller scales only the CHDM 
model is well inside the observational strip. As a more 
quantitative analysis, we estimate, for each model, the 
fraction of observers whose measured bulk fiow has cer- 
tain specific properties. The corresponding probabili- 
ties are also reported in Table |^. In Column 3 we list 
the fraction of observers jTso measuring a bulk flow at 
50/i~^Mpc within la from the BPS value. We choose 
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Table 3. The fractions of observers measuring a particu- 
lar characteristic in different cosmological models. Column 
2: satisfying the Local Group requirements (see text; Tlg)\ 
Column 3: bulk flow within Icr from the Branchini ct al. 
(1995) value at 50 /i^^Mpc (Tt,o)\ Column 4: bulk flow within 
1(T from the Branchini ct al. (1995) value at 150 ft.^-'^Mpc 
(.^l5o); Column 5: bulk flow within Icr from the Laucr & 
Postman (1994) value at ~ lOO/i^^Mpc {Tlp); Column 6: 
bulk flow within 2(t from the Branchini et al. (1995) values 
in the range from 50 to 150/i-lMpc {J^lii)- 



Model 




^50 


.7-150 




7r2(T 
■'all 


SCDM 


0.026 


0.293 


0.365 


0.032 


0.109 


CHDM 


0.053 


0.437 


0.315 


0.019 


0.237 


OCDM02 


< 10-3 


0.017 


< 10-3 


< 10-3 


0.019 


OCDM03 


< 10-3 


0.046 


0.013 


< 10-3 


0.094 


OCDM04 


0.015 


0.383 


0.108 


< 10-3 


0.228 


OCDM06 


0.070 


0.409 


0.300 


0.009 


0.187 


OCDM08 


0.038 


0.343 


0.379 


0.021 


0.151 


ACDM02 


0.052 


0.396 


0.385 


0.047 


0.278 


ACDM03 


0.070 


0.299 


0.403 


0.088 


0.161 


ACDM04 


0.051 


0.274 


0.418 


0.068 


0.141 


ACDM06 


0.038 


0.312 


0.375 


0.042 


0.135 


ACDM08 


0.030 


0.310 


0.382 


0.020 


0.119 



this depth since it is still within the range spanned 
by galaxy catalogues incorporating peculiar velocities. 
Moreover, at 50 ft-^Mpc the BPS bulk velocity over- 
laps with that obtained from the POTENT analysis 
(Dekel 1994) and with that measured using the recent 
da Costa et al. (1995) data. Similarly, Column 4 con- 
tains the probability measured at 150/i-^Mpc (.T^iso). 
Although this is also the limiting depth of the Lauer & 
Postman (1994) sample, we have compared their bulk 
velocity (Vbuik ~ 700 km s-^), which is much higher 
than that of BPS, at a depth of ~ 100/i^^Mpc which 
corresponds to the effective depth of their sample. The 
fraction of observers measuring a bulk flow within la 
from the Lauer & Postman (1994) one is listed in Col- 
umn 5 {Tlp)- Finally we measured the fraction of ob- 
servers whose cumulative bulk velocity is within the 
2(7 strip of the BPS data throughout the range from 
50 to 150/i^^Mpc (Column 6, Tlu). In this case we 
preferred to consider a larger range (2(7) in order to 
have more robust statistics. The comparison with data 
at 50/i-^Mpc and 150/i-^Mpc confirms the visual im- 
pression: the only disfavoured models arc the OCDM 
ones with Oq = 0.2 and 0.3 and, more marginally and 
only at 150 ft-^Mpc, the fio = 0.4 case. Moreover, for 
all models the Lauer & Postman (1994) observation is 
a very low probability event. The same conclusion was 



reached using various techniques (Feldman & Watkins 
1994; Tegmark, Bunn & Hu 1994; JafFe & Kaiser 1994; 
Strauss et al. 1995). Finally, in the open Qq = 0.2 and 
0.3 CDM models, a very tiny fraction of observers mea- 
sures a cumulative bulk velocity in accordance with the 
BPS data in the whole range of scales spanned by data; 
the other models give similar results. 

4.3 Statistics of the velocity distribution 

In order to compare the data properly with the sim- 
ulations, we need to take account of the (fairly large) 
uncertainties in cluster velocity determinations. To this 
end, we convolve simulated cluster velocities using the 
following procedure. Let Vi be the intrinsic line-of-sight 
velocity for the i-th model cluster {i = 1, . . . , 16000) 
and 5vj the velocity error for the j-th randomly selected 
real cluster. The convolved velocity, Vi, of the model 
cluster is obtained by randomly generating it from a 
Gaussian distribution with mean Vi and dispersion 5vj . 
After repeating this operation for all the clusters in a 
simulation, we end up with a sample of error-convolved 
cluster peculiar velocities. 

To estimate the uncertainty in the CVFD for the 
observational data, due to individual velocity errors, we 
resample the real cluster velocities as follows. For each 
cluster having raw velocity v and error 5v, we generate 
a velocity v' which is randomly taken from a Gaussian 
distribution, having mean v and dispersion 5v. A re- 
sampling of the original data set is therefore obtained 
by repeating this procedure for all the clusters. We gen- 
erate 20,000 resamplings of the observational catalogues 
and for each of them we estimate the CVFD. The error 
is estimated as the scatter within this ensemble. The 
CVFD for the composite sample is plotted in Figure 6, 
where for comparison we show also the CVFDs obtained 
considering some of the single samples, namely the Hud- 
son (1994) IRTF sample, the Hudson (1994) _D„-cr sam- 
ple, the Giovanelli (1995) sample and the Branchini et 
al. (1995) reconstructed data (see the following subsec- 
tions for more details about these samples). For sake of 
clarity we plot error bars, corresponding to the 2a scat- 
ter over the observer ensemble, only for the Giovanelli 
(1995) and for the composite samples. 

4..3.1 Results for the composite observational sample 

As a first characterization of the cluster velocity field, 
we compute for our composite sample the rms one- 
dimensional velocity, cr„ = (v^o)^^^- The correspond- 
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ing values for the simulated catalogues are reported in 
Table both before (Column 2) and after (Column 3) 
error-convolution in the manner described above. The 
latter have to be compared with results from the line- 
of-sight observational peculiar velocities, which give 
(T„ — 725 ± 60kms~^. The errors for the convolved 
(Ti, are estimated as the la scatter over an ensemble of 
20,000 samples, each containing 62 randomly selected 
clusters taken from the whole simulation. Note that 
the convolution with observational errors changes sig- 
nificantly the values of av This is particularly true for 
the models with small Qq (and small peculiar motions) 
where a large part of the measured signal comes from the 
error noise. In Column 4 we report the fraction Tcr^ of 
samples measuring values of the rms velocity (t„ within 
the la observational band. As our main result, we note 
that observers in OCDM models with Qq < 0.4 have a 
very low probability of measuring the observed a^ , while 
this value becomes typical if one considers Qq > 0.6. 
Among the ACDM models, only the flo = 0.2 case is 
rather unlikely, while all the other models are in good 
agreement with data. Both the f2o = 1 models are con- 
sistent with the available observations. 

The CVFDs resulting from both real data and sim- 
ulations are compared in Figure 7, where the model 
curves represent the mean over 20,000 samples and the 
observational error bars are 2a uncertainties. Again, 
models with very low density are disfavoured, especially 
for the OCDM case: they produce too few clusters with 
sufficiently high velocity. Moreover, we note that, al- 
though the ACDM08 and SCDM models have the cor- 
rect a-u, they nevertheless appear to have a system- 
atically smaller number of clusters with low velocities 
(< 600kms~^). The other models appear to lie within 
the 2a observational uncertainties for the whole range of 
velocities considered. The discriminatory power of this 
test is strongly reduced by the smallness of the observa- 
tional sample, as well as by the size of the uncertainties 
in the velocity measurements. In fact, as shown in Fig- 
ure 2, this test would be efficient in constraining the 
different models, were it not for the convolution with 
the observational errors which has the efi'ect of partially 
washing out differences between the model CVFDs. 

In order to perform a more quantitative analysis, 
for each model we computed also the fraction of ex- 
tracted samples for which the measured CVFD assumes 
a value within la from the observed data. We are in- 
terested in studying the velocity distribution for values 
corresponding to the low-velocity tail, to the peak and 



to the high-velocity tail of the distribution. Therefore, 
we decided to consider the CVFD at three different val- 
ues for the velocity, i.e. v = 0.5cr„ = 362.5 kms^^, v = 
1(T„ = 725kms~^ and v — 2a^ = 1450 kms~^, where we 
have taken for a^ the observational value; we named the 
corresponding fractions J-q.s, T\ and T2 respectively. 
For our data sample we find P(> 0.5cri,) = 0.537±0.051, 
P(> lcr„) = 0.245±0.041 and P(> 2a-,) = 0.065±0.024. 
The results for the different models, also reported in 
Table confirm our previous qualitative analysis. The 
CVFD at low-velocity rules out the OCDM models with 
< 0.4, and, more marginally, ACDM08 and SCDM 
models; at high velocity the models in disagreement 
with the real data are the open ones with Qq < 0.4 
again, and ACDM02. 

4-3.2 Results for single observational samples 

We consider now three separate samples of cluster pecu- 
liar velocities, in order to study the stability of our pre- 
vious results with respect to the source and the method 
used to obtain the data (i.e. IRTF vs. D„-a). In partic- 
ular we consider: 

(1) The Hudson (1994) Dn-a sample, reported in his 
Table 2. This is a list of 17 objects with distance smaller 
than 8000 km s~^ Unlike the Burstein (1990) sample, 
from which these data originate, Coma is used as the 
calibrating cluster; the reported random error in the cal- 
ibration of the Dn-a relation is 3.7%. As in the previous 
analysis, we prefer to discard the data of Cen45 because 
of the large discrepancy with other determinations. 

(2) The Hudson (1994) IRTF sample, reported in his 
Table 3, a list of 21 objects with a distance smaller than 
8000 km s^^, again originating from Burstein (1990). 
The calibration is made using the nine Aaronson et al. 
(1986) clusters which are located inside the region con- 
sidered. The random error in this calibration is 2%. 

(3) The Giovanelli (1995) sample, an updated ver- 
sion of the Giovanelli et al. (1995) catalogue contain- 
ing 23 cluster peculiar velocities obtained by the IRTF 
method. The calibration comes from those clusters 
which lie at a distance larger than 4000 km s~^. 

The results of the CVFD test for these samples have 
been already shown in Figure 6. In Table |^ we report 
the values of the rms velocity estimated from the line- 
of-sight observational peculiar velocities with the same 
method previously applied to the composite sample. We 
note that, while the results from the IRTF and Dn-a 
Hudson (1994) catalogues are fairly consistent both with 
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Table 4. Results for the composite sample. The different simulation rms line-of-sight velocity values, CTi, (in kms~^), before 
(Column 2) and after (Column 3) convolution with the observational uncertainties of the composite sample. Also reported is 
the fraction Ta^ of the extracted samples (Column 4) having cr^ falling within the la observational range. Columns 5, 6 and 7 
refer to the fractions of extracted samples for which the measured cumulative velocity frequency distribution P{> v) assumes 
a value within la from the observed data, computed for v = O.bay (J-'o.s), v = la^ (jFi) and v = 2cr„ (JF2), respectively. 



Model 


cTy (unconv.) 


(jy (conv.) 




^0.5 


^1 


^2 


Composite (D„-cr & IRTF) 




725 ± 60 










SCDM 


664 


797 ± 78 


0.413 


0.187 


0.133 


0.543 


CHDM 


467 


633 ± 66 


0.293 


0.548 


0.505 


0.276 


OCDM02 


93 


423 ± 61 


< 10"» 


< 10-^ 


0.002 


0.021 


OCDM03 


98 


434 ± 62 


0.001 


< 10-3 


0.002 


0.035 


OCDM04 


197 


465 ± 60 


0.002 


0.013 


0.009 


0.039 


OCDM06 


415 


592 ± 63 


0.123 


0.503 


0.379 


0.183 


OCDM08 


574 


718 ± 72 


0.594 


0.373 


0.390 


0.507 


ACDM02 


247 


492 ± 59 


0.004 


0.066 


0.029 


0.069 


ACDM03 


388 


589 ± 63 


0.115 


0.488 


0.369 


0.173 


ACDM04 


504 


666 ± 68 


0.450 


0.499 


0.523 


0.371 


ACDM06 


570 


720 ± 72 


0.597 


0.358 


0.382 


0.506 


ACDM08 


623 


762 ± 73 


0.546 


0.235 


0.205 


0.558 



each other and with the value obtained from the consol- 
idated catalogue (646 ± 102 kms"^ and 688±82kms-^ 
versus 725 ± 60kms-^), the value of cr„ from the Gio- 
vanelli (1995) sample is much smaller, 356 ± STkms"^. 
The resulting discrepancy has a high confidence level, 
(~ 3a) for both the D„-o- and the IRTF Hudson sub- 
samples. 

In order to compare these results with the cosmo- 
logical models, we convolve the simulated cluster veloci- 
ties with the observational uncertainties for each sample 
separately: the resulting are shown in Table ^. These 
values and their errors are the mean and the la scatter 
over an ensemble of 20,000 samples, respectively. Also 
reported is the fraction Ta^ of those samples which mea- 
sure a value of the rms velocity a^ within the la obser- 
vational band. Once again, the results for the D„-a and 
the IRTF Hudson (1994) samples are very similar and 
in good agreement with the results for the composite 
catalogue: the only models we can reject on the basis of 
this analysis are OCDM with Qo < 0.4 and ACDM02. 

The situation is radically different when the Gio- 
vanelli (1995) data are used. Due to the smaller errors 
reported for this sample, the discriminatory power of the 
test is very high and only two models survive: namely 
OCDM04 and ACDM02. Note that these models are 
both rejected by either the composite or the Hudson 
(1994) samples. 



4-3.3 Results for the BPS reconstructed velocities 

The results of the comparison between models and BPS 
cluster velocities are also shown in Table ^ where we 
report the BPS rms line-of-sight velocity a^ = 522 ± 15 
km , estimated with the same convolution method 
that was previously applied to the observational sam- 
ples. Following the same procedure as before, we find 
that the only models having a non-negligible fraction of 
"good" observers are CHDM, ACDM04 and OCDM06. 
It is interesting to note that now the SCDM model is 
rejected at a high significance level, in agreement with 
the result based on the Giovanelli (1995) data. Note that 
none of the models that satisfy the BPS velocities are 
the same time accepted by the analysis based on the 
Giovanelli data, but are accepted on the basis of the 
Hudson and composite samples. 

It is worth recalling here that the high discrimina- 
tory power of such reconstructed data is due to both 
their large number and their small errors. However, as 
we already mentioned, such errors include observational 
as well as intrinsic uncertainties related to the recon- 
struction procedure but they do not account for possi- 
ble systematic effects related to the validity of the linear 
biasing assumption on which the reconstruction method 
relies. It is wise, therefore, to interpret these results only 
as indications, but indications which are, nevertheless, 
in agreement with the analysis of the composite and 
Hudson data. 
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[table] 

Table 5. Results for single samples. The values of the simulation rms velocities av (in kms~^), evaluated for one-dimensional 
data after convolution with the observational uncertainties, and the fraction J^cry of the extracted samples having cry falling 
within the Icr observational range are reported for the Hudson (1994) Dn-cr catalogue (Columns 2 and 3), for the Hudson 
(1994) IRTF catalogue (Columns 4 and 5), for the Giovanelli (1995) IRTF catalogue (Columns 6 and 7) and for the Branchini 
et al. (1995) reconstructed velocities (Columns 8 and 9). 



Sample 


Hudson D 


'n-cr 


Hudson IRTF 


Giovanelli IRTF 


Branchini et al. 


Objects 


16* 




21 




23 




280 




Model 










Tv 








Observed 


646 ± 102 




688 ± 82 




356 ± 37 




522 ± 15 




SCDM 


786 ± 72 


0.303 


781 ± 73 


0.343 


702 ± 64 


< 10-3 


703 ± 64 


0.003 


CHDM 


625 ± 61 


0.885 


614 ± 60 


0.400 


513 ±50 


0.001 


512 ±50 


0.220 


OCDM02 


407 ± 45 


0.003 


395 ± 47 


< 10-3 


215 ± 27 


0.003 


214 ± 27 


< 10-3 


OCDM03 


418 ± 47 


0.008 


405 ± 50 


< 10-3 


222 ± 27 


0.004 


222 ± 28 


< 10-3 


OCDM04 


454 ± 49 


0.042 


440 ± 49 


< 10-3 


284 ± 29 


0.264 


284 ± 30 


< 10-3 


OCDM06 


581 ± 56 


0.735 


572 ± 56 


0.169 


465 ± 43 


0.008 


464 ± 43 


0.107 


OCDM08 


707 ± 68 


0.727 


700 ± 67 


0.659 


616 ± 57 


< 10-3 


616 ± 57 


0.055 


ACDM02 


482 ± 50 


0.111 


469 ± 50 


0.003 


331 ± 33 


0.672 


330 ± 32 


< 10-3 


ACDM03 


578 ± 55 


0.724 


570 ± 57 


0.163 


458 ± 43 


0.013 


460 ± 44 


0.093 


ACDM04 


654 ± 62 


0.896 


646 ± 62 


0.575 


551 ± 50 


< 10-3 


551 ± 50 


0.199 


ACDM06 


709 ± 67 


0.718 


703 ± 68 


0.650 


613 ± 57 


< 10-3 


612 ±57 


0.059 


ACDM08 


751 ± 68 


0.483 


746 ± 69 


0.513 


668 ± 60 


< 10-3 


667 ± 60 


0.010 



the data of Cen45 is discarded. 



5 DISCUSSION 

This work was devoted to a thorough investigation of 
the velocity field as traced by galaxy clusters. We have 
focusscd our attention on understanding whether such 
large-scale flows can be usefully employed to constrain 
the values of the cosmological parameters fio and Q,a- 

We have used a simulation technique based on the 
Zel'dovich approximation, first testing this method care- 
fully by compaxing its results with a full N-body code. 
We have introduced two different definitions of veloc- 
ity in our simulations. The first, linear, definition (LIN) 
amounts to the assumption that the particle velocities 
are given by the gradient of the linear gravitational 
potential estimated at the initial (Lagrangian) particle 
position. In the second, quasi-linear, definition (ZEL) 
the particle velocity is estimated through the gradi- 
ent of the (Zel'dovich-evolved) potential evaluated at 
the final (Eulerian) particle position. We compared the 
cluster velocity fields corresponding to these two pre- 
scriptions with that obtained from the N-body simu- 
lation by applying different statistical tests, such the 
velocity frequency distribution, the velocity correlation 
function and the bulk flow. Prom this analysis we con- 
clude that the quasi-linear definition (ZEL) for the 
cluster velocities is consistent with N-body results for 



scales 10 — 20 hr^Mpc. This is true even for high- 
normalization spectra, as — 1, when a substantial de- 
gree of non-linearity is present on the cluster mass scale. 

In the framework of the CDM cosmogony, we con- 
sidered models with 0.2 < Ho < 1, both with fiat 

(ACDM) and open (OCDM) geometry. We also con- 
sidered the ^0 = 1 Cold+Hot DM (CHDM) model with 
30% of hot component. For each model, we ran a single 
realization within a large enough box, 960/i-^Mpc on 
a side, to render any effect of cosmic variance entirely 
negligible. 

We first compared predictions of the relevant mod- 
els for both the 3D cluster velocity frequency distribu- 
tion (VFD), P{v), and the velocity correlation function, 
i,v{r). These statistics have been shown to possess an in- 
trinsic discriminatory power between such models. This 
is particularly true for the OCDM models which have 
rather different VFDs and velocity correlations. The dif- 
ferences axe instead smaller in the case of the ACDM 
models, with the only exception of the Q,q — 0.2 case. 
We also found, in agreement with previous analyses (e.g. 
Bahcall et al. 1994a), that all the models display a P{v) 
which is well fitted by a Maxwellian distribution. 

In order to constrain the models we compared their 
predictions with available cluster velocity data. As a 
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first test we computed the probability of having ob- 
servers with the Local Group characteristics (i.e. similar 
local velocity and density contrast). The comparison for 
the bulk flow statistic has been made with the results 
recently obtained by Branchini et al. (1995), based on 
an iterative procedure used to reconstruct the cluster 
real space positions and peculiar velocities. Both anal- 
yses show that the only models which can be ruled out 
at a rather high confidence level are the OCDM ones 
with fio < 0.4. Note that the bulk flow measured by 
Branchini et al. (1995) at 150 /i~^Mpc has a much lower 
amplitude (and much better alignment with the CMB 
dipole) than that obtained by Lauer & Postman (1994), 
which however used cluster peculiar velocities measured 
directly using a distance indicator based on the proper- 
ties of brightest cluster galaxies. We found, in agreement 
with other analyses, that the Lauer & Postman (1994) 
bulk flow is an extremely rare event in most of our mod- 
els, although in the case of ACDM models (particularly 
for those with SIa = 0.3 and 0.4) there is a small but 
non-negligible probability of obtaining the sample am- 
plitude as the Lauer & Postman bulk flow. 

As a further test, we computed the ID cumulative 
velocity frequency distribution (CVFD) and the veloc- 
ity dispersion Gv We constructed a composite catalogue 
by collecting observed cluster velocities from different 
data sets, based on both Dn-a and IRTF methods (see 
Table ^). The flnal observational sample used in the 
analysis contains 62 objects. We took into account the 
large observational uncertainties by convolving the sim- 
ulated cluster velocities with the error reported in the 
source catalogues. This convolution broadens the veloc- 
ity distributions and increases the values of , thus re- 
ducing the differences between the models. In this way 
the possibility of constraining the Q,o and f^A param- 
eters is largely reduced. The low-57o models are again 
disfavoured, particularly in the case of open models, for 
which f2o > 0.4 is required to fit the data; the ACDM02 
model is also in trouble, while ACDM08 and SCDM 
have systematically fewer clusters with high velocity 
than in the observational sample. 

In order to check whether our composite sample 
is affected by biases associated with its heterogenous 
composition, we considered also three smaller samples 
directly taken from the literature: the Dn-a and the 
IRTF Hudson (1994) catalogues and the IRTF Gio- 
vanelli (1995) catalogue. The paucity of the data and 
the large uncertainties affecting the two Hudson sam- 
ples yield less stringent results, but in agreement with 



the analysis of the composite catalogue (cf. Tables |j 
and Note that the two Hudson samples contain data 
obtained from different distance indicators {D„~cr and 
IRTF relations). The similarity of the results therefore 
seems to indicate that there is not an intrinsic bias 
introduced by the technique adopted in obtaining the 
data. On the other hand, the results from the Giovanelli 
(1995) sample are completely different: this catalogue is 
dominated by clusters with very small velocities and the 
resulting velocity dispersion is consequently also very 
small. In this case, the only models in agreement with 
the data are OCDM04 and ACDM02, which were ex- 
cluded by the analysis of the composite and of the Hud- 
son samples. 

The discrepancy between the Giovanelli data and 
all other studies was already noted by the authors them- 
selves (see, e.g., Giovanelli et al. 1995). They claim that 
their results are more reliable for the following reasons: 

(1) they used a homogenous galaxy sample, which pro- 
vides a full coverage of the sky; 

(2) they adopted a different TuUy-Fisher template re- 
lation, obtained from an extensive study of clusters; 

(3) they applied an internal extinction correction 
which depends on the luminosity and allows for larger 
flux corrections. 

On the other hand, a possible bias for the Gio- 
vanelli sample might have arisen from their estimate of 
the IRTF relation zero-point. They adopted that of the 
mean relation obtained from all clusters with cz > 4000 
km s~^ in the CMB frame, since, as they argue, at 
such distances small peculiar velocities translate into 
small magnitude offsets, and that possible bulk flows do 
not affect the mean IRTF relation since the clusters are 
isotropically distributed. However, since the number of 
available clusters is small (~ 15), even if they are ran- 
dom samplers of an isotropic distribution, they could in- 
troduce shot -noise deviations from the 'isotropic' zero- 
point value. A second concern is related to the use of the 
CMB rest-frame. If there is a bulk flow, in which the LG 
participates, then, by using this frame the sampled vol- 
ume becomes elongated along the direction of the bulk 
flow, compared to a spherically symmetric distribution. 

We have also analysed the CVFD for the BPS re- 
constructed cluster velocities, which represents a com- 
pletely independent data set. It turns out that such ve- 
locities pose the most stringent constraints on the dif- 
ferent cosmological models (cf. Table Since, however, 
the BPS velocities are not truly observational but de- 
pend on various assumptions we choose not to present 
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the above results as final and firm conclusions. In par- 
ticular, the assumption of a linear biasing relationship 
which is independent of spatial position, may well not 
be valid. Our ignorance about the possible cluster-by- 
cluster variation of the biasing parameter would increase 
the individual reconstructed velocity uncertainties. For 
these reasons, we prefer at present to consider these re- 
sults only as circumstantial evidence in support of infer- 
ences made on the basis of the composite and Hudson 
data. 



6 CONCLUSIONS 

The first definite conclusion of this work is that our 
method to simulate the cluster velocity field, which is 
based on an optimized version of the Zel'dovich approx- 
imation, is easy to implement, cheap from the computa- 
tional point of view and fully reliable, as shown by the 
comparison with the results of a full N-body code. 

As for the comparison between observational data 
and numerical simulations on cluster peculiar velocities 
is concerned, our main findings can be sketched as fol- 
lows. 

(a) The bulk fiow analysis performed on the Branchini 
et al. (1995) data and the probability of reproducing the 
Local Group characteristics rule out only open CDM 
models with SIq < 0.4. We also confirm that the Lauer & 
Postman (1994) high bulk velocity, at an effective depth 
of ~ 100/i~^Mpc, is an extremely unlikely event in all 
the r^A ~ models. However, in Qa models there 
is a small but non-negligible probability of reproducing 
this observation. 

(b) The VFD analysis based on a composite sample of 
62 cluster velocities, obtained both with IRTF and Dn~ 
a methods, rules out open CDM models with f2o < 0.4 
and the ACDM model with Slo = 0-2. Standard CDM 
and ACDM with f2o = 0.8 are disfavoured, although 
only marginally. Similar results are found by applying 
the same test on the two Hudson (1994) catalogues, 
based on the IRTF and D„-a relations, respectively: 
no systematic bias appears to be introduced by the dif- 
ferent distance indicators. 

(c) The application of the VFD test to the IRTF Gio- 
vanelli (1995) sample gives completely different results. 
The smaller velocities and errors reported in this sam- 
ple are such that the only surviving models are the open 
CDM with Qo = 0.4 and the ACDM with = 0.2. 

(d) We found also that the above two models, which 
are consistent with the Giovanelli data, are strongly dis- 



favoured by an analysis of the Branchini et al. (1995) 
reconstructed velocities, in agreement with the interpre- 
tation of the Hudson and composite sample analyses. 

As a concluding remark, we would like to stress the 
importance of identifying the reasons for the puzzling 
discrepancy between different data sets. Our analysis 
has shown that a small, but accurate, sample of cluster 
peculiar velocities, like that of Giovanelli (1995), can put 
rather stringent constraints on different DM models. We 
therefore consider the analysis presented here as proving 
the potential of cluster velocity considerations, and look 
forward to the availability of newer and more precise 
cluster velocity measurements. 
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FIGURE CAPTIONS 

Figure 1. Comparison of tiie projected velocity fields 
traced by clusters in a N-body simulation (PM, up- 
per row) and in Zel'dovich simulations with two dif- 
ferent definitions for the velocity: ZEL (central row) 
and LIN (bottom row). The simulations have a box- 
size of 320/i"^Mpc and the same CHDM initial condi- 
tions. The slice is 60/i~^Mpc thick. A smoothing by a 
Gaussian filter with width 20/i~^Mpc is applied. Two 
different stages are shown: as = 0.67 (left column) and 
(78 = 1 (right column). 

Figure 2. Comparison of the results of the application 
of different statistical tests to the cluster velocity fields 
obtained in a N-body simulation (PM, solid lines) and 
in Zel'dovich simulations with two different definitions 
for the velocity: ZEL (dotted lines) and LIN (dashed 
lines). Upper row: the 3D velocity frequency distribution 
P{v). Central row: the velocity correlation function f„ 
(in km s^^) as a function of the distance r (in ft^^Mpc). 
Bottom row: the bulk flow velocity Vbuik (in km s"'^) as 
a function of the distance r (in /i~^Mpc). Two different 
stages are shown: era = 0.67 (left column) and cts = 1 
(right column). 

Figure 3. The frequency distribution P{v) for the 3D 
cluster velocity field for different models: SCDM and 
CHDM models (top left), OCDM models (top right), 

ACDM models (bottom left). The dashed lines in the 
top left panel show the Mcixwellian curves fitting the 
distributions. 

Figure 4. The velocity correlation function ^„ (in km 

s~^) as a function of the distance r (in h~^Mpc) for 
different models: SCDM and CHDM models (top left), 
OCDM models (top right), ACDM models (bottom 
left). 

Figure 5. The bulk fiow velocity Vbuik (in km s~^) as 
a function of the distance r (in h~^Mpc) for different 
models: SCDM and CHDM models (top left), OCDM 
models (top right), ACDM models (bottom left). The 
dashed region represents the la strip allowed by the 
Branchini et al. (1995) data. 

Figure 6. The ID cumulative velocity frequency distri- 
bution P(> v) of our composite catalogue (solid line) 
is compared to that of the Hudson (1994) IRTF sample 
(short-dashed line), of the Hudson (1994) Dn a sam- 
ple (long-dashed line), of the Giovanelli (1995) sample 
(dotted line) and of the Branchini et al. (1995) recon- 
structed data (dotted-dashed line). Error bars, plotted 
only for the composite and Giovanelli samples, corre- 
spond to the 2cr scatter over 20,000 resamplings. 



Figure 7. The ID cumulative velocity frequency distri- 
bution P(> v) for different models compared with that 
resulting from our composite observational catalogue 
(filled circles with 2(7 uncertainties obtained as scatter 
over an ensemble consisting of 20,000 resamplings of the 
original data). The model curves are the mean of 20,000 
samples of 62 randomly selected clusters whose veloc- 
ities are convolved with the observational errors. The 
panels refer to SCDM and CHDM models (top left), 
OCDM models (top right) and ACDM models (bottom 
left). 
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